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In this article I review the main theoretical problems that are posed by the highest energy end of the observed 
cosmic ray spectrum, stressing the importance of establishing their composition in order to decide between pro- 
posed scenarios. I then discuss the possibilities that are opened by the detection of inclined showers with extensive 
air shower arrays. Recent progress in modelling magnetic deviations for these showers has allowed the analysis 
of inclined showers that were detected by the Haverah Park experiment. This analysis disfavours models that 
predict a large proportion of photons in the highest energy cosmic rays and open up new possibilities for future 
shower array detectors particularly those, like the Pierre Auger Observatory, using water Cerenkov detectors. 



1. HIGH ENERGY COSMIC RAYS 

Cosmic rays are elementary particles arriving 
at the Earth from outside that were discovered in 
the beginning of the 20th century as one of the 
main sources of natural radiation. The cosmic ray 
spectrum has been observed as a continuum at all 
energies since their discovery. Throughout this 
period cosmic rays have always been the source 
of the highest energy elementary particles known 
to mankind, and for this reason they have given 
birth to particle physics. The high energy tail of 
the spectrum as it is known today corresponds 
to energies up to 3 10^° eV and rates of a few 
particles per km^ per century. 

It is remarkable that the cosmic rays have 
a quite featureless power law energy spectrum 
which decreases as approximately the cube of the 
primary energy. For energies above the few hun- 
dred TeV the observed flux necessarily requires 
techniques that take advantage of the extensive 
air showers that the arriving particles develop as 
successive secondary particles cascade down into 
the atmosphere. Shower measurements allow the 
reconstruction of the arrival directions and the 
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shower energy but the nature of the primary par- 
ticle is extracted by a number of indirect meth- 
ods. 

For energies above few tens of GeV the detected 
particles, mainly protons, have arrival directions 
with a remarkably isotropic distribution. This is 
understood in terms of diffusive propagation in 
the galactic magnetic fields. As the energy rises 
above a given value that depends on the charge 
of the particle, propagation in the Galaxy should 
cease to be diffusive. Such high energy particles 
are expected to be extragalactic. 

The observation of high energy cosmic rays 
has been recently reviewed by Nagano and Wat- 
son who have shown that there is very good 
agreement between different experiments includ- 
ing the low and high energy regions of the spec- 
trum. There is increasing evidence for a different 
component of the high energy end of the cosmic 
ray spectrum |^. Combining data of five dif- 
ferent experiments, AGASA, AKENO, Haverah 
Park, Stereo Fly's Eye and Yakutsk, Nagano and 
Watson conclude that there is a clear signal of 
a change of the spectral slope in the region just 
above 10^® eV 0. Composition studies have also 
given indications that there is a change to light el- 
ement composition for energies above ~ 10^^ eV 
S although this conclusion is model dependent 
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to some extent g]. Also the small anisotropy ( 
4%) of 10^^ eV cosmic rays in the direction of the 
galactic anticenter detected with AGASA disap- 
pears at higher energies Q. 

The highest energy events detected present a 
serious challenge to theory and little is known 
about their origin. If they are protons they should 
attenuate in the Cosmic Microwave Background 
(CMB) over distances of order 50 Mpc. Such at- 
tenuation was predicted to appear in the cosmic 
ray spectrum as a cutoff, the Greisen-Zatsepin- 
Kuz'min (GZK) cutoff, just above 4 10^° GeV 
Q . If they are photons or iron nuclei it turns out 
that interactions with the radio and the infrared 
backgrounds are respectively responsible for at- 
tenuations over similar or even shorter distances. 
No such features are seen in the observed cosmic 
ray spectrum. If they are produced sufficiently 
close to us to avoid the cutoff then the arriving 
particles should be pointing to their sources. This 
seems difficult to accommodate because there are 
very few known astrophysical sources capable of 
reaching the observed energies and on the other 
hand there is little evidence for the anisotropy 
that would result. 

This article firstly discuses the problem pre- 
sented by the high energy end of the cosmic ray 
spectrum with emphasis in the role of composi- 
tion. Then it outlines new progress made in un- 
derstanding different features of inclined showers 
illustrating how these showers can contribute to 
the composition issue reviewing the results ob- 
tained by a recent analysis of the inclined data in 
Haverah Park. 

2. ORIGIN: AN UNSETTLED ISSUE 

The discovery of events with energies above 
10^0 eV (100 EeV) dates back to the 1960's, to 
the early days of air shower detection experiments 
1^ . Since then they have been slowly but steadily 
detected by different experiments as illustrated 
in Fig. |l|. Now there is little doubt about the 
non observation of a GZK cutoff, with over 17 
published events above 10^° eV and five prelim- 
inary new events from HiRes 0. On the con- 
trary the data suggests that the spectrum contin- 
ues smoothly within the statistical errors, possi- 



bly with a change of slope. On the other hand the 
data show no firm evidence of anisotropy but the 
significance of such studies is even more limited 
by the poor statistics. 
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Figure 1. Accumulated events of energy exceed- 
ing 10^*^ eV plotted function of time as de- 
tected by different experiments: Volcano Ranch 
(VR), Haverah Park (HP), Horizontal Air Show- 
ers in Haverah Park (HAS-HP), Yakutsk, Fly's 
Eye and AGASA. 



Both the details of the spectrum at the cut- 
off region and the extent to which the arrival di- 
rections of these particles cluster in the direction 
of their sources are very dependent on a num- 
ber of unestabhshed issues. These include the 
source distribution, the distance of the nearest 
sources, their emission spectra, the intervening 
magnetic fields and of course on the nature of 
the the cosmic rays themselves or composition. 
If these particles are nuclei or photons the ob- 
servational evidence is suggesting that these are 
coming from relatively nearby sources compared 
to the 50 Mpc scale. The conclusive power of 
observations is however strongly limited by both 
the poor statistics and a complex interrelation of 
hypotheses, but the situation is bound to change 
in the immediate future with a new generation of 
large aperture experiments, some like HiRes M 



already in operation, others in construction 
and many others in planning |]9|,p^ . The complex 
puzzle that connects particle physics, magnetic 
fields, and cosmic rays has attracted the atten- 
tion of many fields in physics. 

In a conventional approach these particles 
would be nuclei as the bulk of the cosmic ray 
spectrum which are accelerated through stochas- 
tic acceleration as suggested by Fermi in 1949. 
This happens every time charged particles cross 
interfaces between regions that have astrophysi- 
cal plasmas with different bulk motions, such as 
shock fronts. Transport is assumed to be diffusive 
in the plasma's magnetic field and on average in 
these processes a very small fraction of the bulk 
plasma kinetic energy is transferred as a boost 
to the individual particles, that typically end up 
with a power like spectrum. Acceleration of a par- 
ticle of charge Ze to an energy E is strictly limited 
by dimensional arguments to objects that are suf- 
ficiently large or have sufficiently large magnetic 
fields. Basically for a particle with momentum p 
to be able to undergo such a boost, propagation 
must be diffusive, or equivalently the accelerator 
region L must be larger than the Larmor radius 
of the particle, i?, in its characteristic magnetic 
field B: 

R = — ^ < L E~pc< ZeBcL (1) 
ZeB 

The requirement is well known by accelerator de- 
signers and is the ultimate reason for their high 
cost. It turns out that few of the known as- 
trophysical objects satisfy the minimum require- 
ments to accelerate particles to 10^" eV. This is 
conveniently illustrated in a plot first conceived 
by Michael Hillas ||T^ which is reproduced in 
Fig. ^. A number of possible scenarios are be- 
ing discussed; they imply acceleration in some ob- 
jects including young pulsars. Gamma Ray Bursts 
(GRB), our own galaxy, active galaxies and the 
local group of galaxies [|2j. The power sup- 
ply needed to keep the observed cosmic rays at 
the highest energies is consistent with the known 
power and distributions of these objects ||^. 

It is difficult to explain the observed flux spec- 
trum in this conventional approach. A solution 
in which particles are accelerated nearby has dif- 
ficulties because there are very few objects which 
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Figure 2. Hillas Plot of the typical size of a 
possible accelerator L versus its magnetic field 
strength B. From upper left to lower right 
the astrophysical objects correspond to Neutron 
Stars, White Dwarfs, Sun Spots, Magnetic Stars, 
Gamma Ray Bursts, Active Galactic Nuclei, Inter 
Planetary Space, Supernovae Remnants, Radio- 
galaxy Lobes, Galactic Disk, Galactic Halo, Clus- 
ters of Galaxies and the Inter Galactic Medium. 
Also shown is the point corresponding to the 
largest accelerator in planning LHC. The straight 
lines represent the limits given by Eq (13) for pro- 
tons (full). Iron nuclei (dashed) and for protons 
assuming a 10% efficiency (dots). 



are capable of accelerating particles to the maxi- 
mum observed energies. Moreover many such ob- 
jects are either too large or too distant for the 
cosmic ray spectrum detected at the Earth not 
to show the predicted GZK cutoff. If the sources 
were to be galactic no absorption cutoff would be 
expected but some spectral features are predicted 
for primary protons that are produced at a dis- 
tance of more than a few Mpc. On the other hand 
the non observation of anisotropy complicates the 
puzzle, because the location of the possible accel- 
erators in our vicinity is pretty well known. Pri- 
mary protons having energies in the 10^'' eV range 
are expected to be little deviated in the galactic 
magnetic fields. Our knowledge of extragalactic 
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magnetic fields is poor but bounds on extragalac- 
tic magnetic fields also imply that the deviations 
of protons produced in the few Mpc range are not 
large. There are however possible configurations 
of the extragalactic magnetic fields that could ex- 
plain many of the ultrahigh energy events as com- 
ing from a single source jl^ . The issue is far from 
being resolved and knowledge about composition 
is bound to play a crucial role for future progress 
in understanding. 

Motivated by particle physics beyond the stan- 
dard model, many alternatives have been pro- 
posed that avoid acceleration and others that pos- 
tulate different particles or different interactions. 
These include annihilation of topological defects 
created in the early universe, heavy relics that 
survive from the primeval bath, non thermal par- 
ticles that couple to gravity, or wimpzillas and an- 
nihilation of relic neutrinos with messenger neu- 
trinos coming from remote places . As regards 
composition two large categories of possible sce- 
narios can be made namely those in which the 
observed particles are accelerated and those in 
which they are decay products of other particles. 
These two classes differ greatly in composition. 
A knowledge of composition is doubly important 
because firstly it may decide between these two 
classes of solutions and secondly because it would 
simplify the task of interpreting anisotropy mea- 
surements. 

3. COMPOSITION 

The models that depend on acceleration can 
reach higher energies if the accelerated particles 
have large charge Z. This shows as a different re- 
striction line in Fig ^. The relative composition 
of different nuclei resulting from such a scenario 
will depend on the local abundances of the differ- 
ent nuclei and on the energy. Depending on dis- 
tance to the source and the surrounding environ- 
ment there may be energy losses, absorption and 
the production of secondary particle fluxes. For 
instance in Active Galactic Nuclei (AGN) mod- 
els the accelerated protons are expected to inter- 
act with ambient light or matter to produce pi- 
ous that decay into photons and neutrinos. The 
neutrinos can reach the Earth unattenuated and 



provide a signature of proton acceleration. Un- 
less the environment becomes opaque to protons 
the relative fluxes of neutrinos and protons that 
reach the Earth should be a number of order one 
or smaller, just because neutrinos are secondaries 
with repect to protons. The relative fluxes of pho- 
tons and protons would be similar to neutrinos or 
smaller depending on the photon absorption both 
at the source and during transport to Earth. Ra- 
tios of the same order of magnitude would apply 
to most acceleration models. 

Most of the non accelerating alternatives pos- 
tulate the cosmic rays are products of the de- 
cay of other more massive particles produced by 
different mechanisms. Typically these particles 
of mass of 1024-25 eV (often an X particle) de- 
cay into standard model particles which eventu- 
ally fragment into hadrons, mostly pions and a 
small fraction of order 3% of nucleons. While neu- 
tral pions decay into photons charged pions decay 
into neutrinos. Fragmentation processes, known 
form accelerator experiments and extrapolated to 
the high energies, become the common reference 
point for these mechanisms. For this reason all 
these models share a very similar composition 
dominated by photons and neutrinos which typ- 
ically are about ten times more numerous than 
nucleons at the production site. 

Depending on the source distribution the rela- 
tive fluxes of these particles are modified through 
their interactions with the background radiation 
fields. The neutrinos are the particles that pre- 
serve their production spectrum without being at- 
tenuated. Protons get attenuated in few tens of 
Mpc in the cosmic microwave background, (the 
GZK cutoff), while photons are attenuated al- 
ready in few Mps mainly through pair production 
in the radio background. As a result the ratio 
of neutrinos to protons can in principle become 
higher at the Earth than when they are produced 
if the sources are quite distant or cosmologically 
distributed. Many of the proposed mechanisms 
are expected to cluster in our galactic halo. This 
possibility is receiving a lot of attention because it 
would provide a relatively natural explanation for 
the absence of the GZK cutoff. In that case how- 
ever the sources will be quite near and the ratio of 
photons to nucleons should be expected to be of 
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order 10, close to its value at production. Other 
sources are not expected to cluster and hence the 
photon to nucleon ratio is expected to drop to 
values close to one. The ratio of photons to nu- 
cleons depends on the source distribution and is 
rather sensitive to clustering. 

4. INCLINED SHOWER FEATURES 

Most air shower detectors in existence consist 
on arrays of particle detectors that sample the ex- 
tensive air shower front as it reaches the ground. 
Multiple particle production takes place in the 
successive high energy interactions produced as 
the shower penetrates the medium. As a result 
the number of particles in the shower front in- 
creases exponentially. When the average particle 
energy in the front becomes too low for multiple 
particle production the shower reaches it maxi- 
mum number of particles. The development of 
these showers is typically governed by the radi- 
ation length in the material which is of order 
36 g cm^'^ in air and shower maximum, which 
is only logarithmically dependent on the primary 
particle energy, occurs at a couple of thousand 
meters for vertical showers of energies of order 
1020 eV. 

Vertical showers are close to shower maximum 
when reaching the Earth's surface, have pretty 
good circular symmetry and are less affected by 
the Earth's magnetic field. It is thus not sur- 
prising that air showers have traditionally been 
studied at close to vertical incidence, typically for 
zenith angles below 45°, in summary because it 
is much simpler. Moreover in most extensive air 
shower arrays the particle detectors are oriented 
to have maximum collection area for vertical inci- 
dence. Since these detectors are often scintillator 
sheets, they tend to become very inefficient for 
very inclined showers. 

As the zenith angle increases the traversed at- 
mospheric depth rises from 1000 to close to 36000 
g cm~2. As a result the shower maximum is 
reached in the upper layers of the atmosphere and 
most of the shower is absorbed before reaching the 
ground. It has been known for a long time that 
weakly interacting particles such as neutrinos can 
induce close to horizontal air showers deep in 



the atmosphere with particle distributions that 
are quite similar to vertical showers ||lj,|l^. Air 
shower array detectors looking in the close to hor- 
izontal direction can thus be sensitive to high en- 
ergy neutrino fluxes |l6| . In fact most bounds on 
neutrino fluxes have already been obtained from 
air shower experiments [p7| , p^ . 

The original motivation of studying inclined 
showers was to understand the cosmic ray back- 
ground to the neutrino induced showers. Al- 
though the electromagnetic part of the air shower 
induced by an inclined cosmic ray is indeed ab- 
sorbed before reaching ground level, the shower 
front however also contains muons which are 
mainly produced by charge pion decay when the 
primary particle is a hadron. These muons do 
travel practically unattenuated all the slant at- 
mospheric depth and produce density patterns on 
the ground that are much affected by the Earth's 
magnetic field. It has recently become quite clear 
that such inclined showers can be analysed. This 
not only nearly doubles the aperture of any air 
shower array but, when combined with vertical 
measurements, it has a remarkable potential for 
the study of primary composition jl^ . 

Much development in this field has been possi- 
ble by the modelling of the muon density patterns 
produced by inclined showers under the influence 
of the Earth's magnetic field |^ . The lateral dis- 
tributions of muons in inclined showers can be un- 
derstood in terms of a simple model | |2l[ | in which 
the magnetic field is firstly neglected. The model 
stresses two important facts that have been ex- 
tensively checked with simulations in the absence 
of a magnetic field |2^: Most of the muons in 
an inclined shower are produced in a well defined 
region of shower development which is quite dis- 
tant from the ground and the lateral deviation of 
a muon is inversely correlated with its energy. 

Indeed most of the fundamental properties of 
these inclined showers are governed by the dis- 
tance and depth travelled by the muons. It is 
remarkable that the average slant distance trav- 
elled by the muons is of order 4 km for vertical 
showers, becomes 16 km at 60° and continues to 
rise as the zenith angle rises to reach 300 km for a 
completely horizontal shower. This distance plays 
a crucial role as a low energy smooth cutoff for 
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the muon energy distribution. For inclined show- 
ers the muons must have much more energy at 
production to reach ground level without decay- 
ing than in the vertical case. Both the travel time 
and the muon energy loss become relevant. 

The model simply assumes that all muons are 
produced at a given altitude d with a fixed trans- 
verse momentum pj_ that is uniquely responsible 
for the muon deviation from shower axis. In the 
transverse plane to the shower at ground level the 
muon deviation, f, is inversely related to muon 
momentum p. The density pattern has full cir- 
cular symmetry when there is no magnetic field. 
When the magnetic field effects are considered 
the muons deviate a further distance Sx in the 
perpendicular direction to the magnetic field pro- 
jected onto the transverse plane B±, given by: 

S.^'-^^'-l^r- = af, (2) 

2p p±_ 

where in the last equation _Bx is to be expressed 
in Tesla, d in m and p± in GeV. As the muon devi- 
ations are small compared to d they can be added 
as vectors in the transverse plane and the muon 
density pattern is a relatively simple transform 
of the circularly symmetry pattern. The muon 
patterns in the transverse plane can be projected 
onto the ground plane to compare with data as 
well as standard simulation programs. 

Eq. 1^ is telling us that all positive (negative) 
muons that in the absence of a magnetic field 
would fall in a circle of radius f around shower 
axis, are translated a distance 5x to the right 
(left) of the i?_L direction. The dimensionless 
parameter a measures the relative effect of the 
translation. For small zenith angles d is relatively 
small and a << 1 so that the magnetic effects are 
also small, and results into slight elliptical shape 
of the isodensity curves. 

For high zeniths however a > 1 the magnetic 
translation exceeds the deviation the muons have 
due to their p_L. In this case shadow regions with 
no muons are expected in the muon density pro- 
files. For an approximate p±^ ^ 200 MeV and 
B± = 40 this happens when d exceeds a dis- 
tance of order 30 km, that is for zeniths above 
^ 70°. These shadow regions in the transverse 
plane are indeed an outstanding feature of the 
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Table 1 

Relationship between muon number and primary 
energy for proton and irons in two hadronic mod- 
els (see equation H). 



ground density profiles at high zeniths as seen in 
the simulations. 

The simple model can be actually generalized 
to account for muon energy distributions as a 
function of distance to shower axis, and improved 
using the correlation between the average muon 
energy and the distance to shower axis as ob- 
tained in dedicated simulations. When all this 
is done the obtained muon density patterns are 
shown to be accurately refiect those obtained with 
simulations and this proves to be a very useful 
tool for the study of inclined showers. 

For each zenith angle the primary particle en- 
ergy sets the normalization of the particle den- 
sities. For proton primaries the total number of 
muons in the shower scales with the proton en- 
ergy E as: 

N = Nref (3) 

where /? is a constant. It is remarkable that the 
shape of the lateral distribution of the muons 
does not significantly change for showers of energy 
spanning over three orders of magnitude. The 
same happens for heavier nuclei with slightly dif- 
ferent parameters. The results are slightly model 
dependent. Two alternative hadronic interaction 
models have been compared, the Quark Gluon 
String Model (QGSM) and SIBYLL to give also 
the same behaviour with also different parame- 
ters. Table |l| illustrates these effects. 

As a final result the muon distributions can 
be represented by continuous functions which are 
analytically obtained once we know the main fea- 
tures of a shower in the absence of magnetic field. 
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In practice this implies that only different zenith 
angles have to be simulated. Different azimuths 
are obtained by adequate transformations of the 
showers without magnetic deflections. The algo- 
rithm is fast and allows detector simulation and 
also event by event reconstruction of data ob- 
tained by air shower experiments. 

5. PHOTON COMPOSITION 

This powerful technique has been used to anal- 
yse the inclined shower data obtained in the Hav- 
erah Park array. The Haverah Park detector was 
a 12 km^ air shower array using 1.2 m deep wa- 
ter Cerenkov tanks that was running from 1974 
until 1987 in Northern England which has been 
described elsewhere [Q. It is possibly the most 
appropriate detector for this study because the 
water Cerenkov tanks have a uniquely large cross 
section to sample shower fronts of horizontal air 
showers. Moreover the Cerenkov technique gives 
larger signals for muons than for electrons simply 
because the muons have typically larger energies 
and travel through the whole detector. 

A careful study has been made of the energy 
deposition of signal in water Cerenkov tanks by 
horizontal muons using conventional simulation 
programs for this purpose p^ . A number of ef- 
fects have to be considered to interpret the ob- 
served data. Inclined particles can produce light 
that falls directly into the phototubes without be- 
ing reflected in the tank walls. Horizontal muons 
produce more signal through delta rays because 
on average they have higher energies than in ver- 
tical showers. There is a significant signal de- 
posited by electromagnetic particles that arise 
mainly through muon decay. Finally the higher 
energy muons are more likely to deposit more en- 
ergy in the tanks because of catastrophic energy 
losses. The event rate as a function of zenith an- 
gle has been simulated with careful treatment of 
all these effects using the muon distributions ob- 
tained as described in the previous section. The 
qualitative behaviour of the registered rate is well 
described in the simulation and the normalization 
is also shown to agree with data to better than 
30% using the measured cosmic ray spectrum for 
vertical incidence, assuming proton primaries and 
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Figure 3. Density maps of two events in the plane 
perpendicular to the shower axis. Recorded muon 
densities are shown as circles with radius propor- 
tional to the logarithm of the density. The de- 
tector areas are indicated by shading; the area 
increases from white to black as 1, 2.3, 9, 13, 34 
m^. The position of the best- fit core is indicated 
by a star. Selected densities are also marked. The 
y-axis is aligned with the component of the mag- 
netic field perpendicular to the shower axis. 



using the QGSM model |§. 

More impressive are the results of fits of the 
models for muon densities to the observed parti- 
cle densities sampled by the different detectors on 
an event by event basis. The nearly 10,000 events 
recorded with zenith angles above 60° have been 
analysed for arrival directions, impact point and 
primary energy in the assumption the primaries 
are protons. A complex sequence of arrival direc- 
tion and density fits is performed to minimize the 
effect of correlations between energy and arrival 
directions. 

The analysed date is subject to a set of quality 
cuts: the shower is contained in the detector (dis- 
tance to core less than 2 km), the probability 
of the event is greater than 1% and the downward 
error in the reconstructed energy is less than 50%. 
These cuts ensure that the events are correctly re- 
constructed and exclude all events detected above 
80° . Examples of reconstructed events compared 
to predictions are illustrated in Fig. ^. Two new 
events with energy exceeding 10^" eV have been 
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revealed. The results have been compared to a 
simulation that reproduces the same fitting pro- 
cedure and cuts using the cosmic ray spectrum 
deduced from vertical air shower measurements in 
reference ^ . The agreement between the integral 
rate above 10^^ eV measured and that obtained 
with simulation is striking when the QGSJET 
model is used. Sibyll leads to a slight underes- 
timate [|9|. 

The universality of the muon lateral distribu- 
tion function is very powerful and once the equiv- 
alent proton energy is determined for all events, 
the corresponding energies in the assumption that 
the primaries are iron nuclei (photons) can be ob- 
tained multiplying the proton energy by a factor 
which is ~ 0.7 (6) for 10^^ eV. As a result when 
a photon primary spectrum is assumed the simu- 
lated rate seriously underestimates the observed 
data by a factor between 10 and 20. A fairly 
robust bound on the photon composition at ul- 
tra high energies can be established assuming a 
two component proton photon scenario. The pho- 
ton component of the integral spectrum above 
1019 eV (4 10^9 eV) must be less than 41% (65%) 
at the 95% confidence level. Details of the anal- 
ysis are presented in |l^ . 

The results of this method when applied to a 
first analysis of inclined showers produced by cos- 
mic rays above 10^^ eV demonstrates that the 
study of inclined showers not only can double the 
acceptance of air shower arrays but it can be a 
very useful tool for the study of photon composi- 
tion. 
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